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Abstract. It is well known that diffuse 
gamma-rays are produced by collisions between 
Galactic cosmic rays and interstellar matter 
(Hayakawa Hutchinson Morrison hypothesis). We 
estimate the contribution of high energy protons 
produced by stellar flares to the diffuse component 
of Galactic gamma-rays. Each star in the Galaxy 
regularly produces flares that generate high energy 
protons in the energy range from 100 MeV to 1 TeV. 
These protons are a likely source of cosmic rays. 
We estimate that the Galactic bulge can constantly 
generate 2.7 x 10 14 A of high energy proton current. 
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I. Introduction 

According to the Hayakawa— Hutchinson— Morrison 
hypothesis, the diffuse component of Galactic gamma- 
rays is produced by the collision of high energy cosmic- 
rays with interstellar matter [1-3], while according to 
long-held theories, super novae [4] or the interstellar 
medium [5] are possible sources of cosmic rays [6] 
(See Figre 1). In this paper, we show the other possible 
source of cosmic rays; the diffuse component of Galactic 
gamma-rays originate from the flare particles of stars 
throughout the Galaxy. 

As an analogy for gamma-rays, we consider the model 
of solar wind. According to current theory, solar wind is 
formed by the cumulative effect of micro and nano flares 
on the solar surface [7-8]. The duration of these micro 
and nano flares is short and their scale is small. However, 
these small solar flares occur very frequently and hence 
a continuous high speed solar wind is maintained with 
a constant intensity [9-10]. 

With reference to this model, we propose that the 
origin of the diffuse component of Galactic gamma- 
rays may be attributed to continuously produced high 
energy protons originating from stellar flares. In the 
Galactic bulge, there are about 1.6 x 10 10 stars. (We call 
this number flA ). Therefore, although stellar flares are 
small, they occur frequently in each star, and hence their 
emissions create a steady high energy proton current 
with constant intensity. Thus, the bulge of our Galaxy 
can be likened to a powerful cosmogenic DC current 
generator. 

Occasionally a star in the bulge releases a large 
flare, emitting high energy particles beyond Ep >100 
GeV. These high energy particles collide with the stellar 
atmosphere or interstellar matter and produce gamma- 
rays with energy E 7 >10 GeV. In the next section, we 



show numerically that the present model can explain the 
intensity of the diffuse component of Galactic gamma- 
rays. 

II. Fundamental Data used for Calculations 

Before describing our model numerically, we show 
here several numerical values that are necessary for 
defining the model. 

(1) It is well known that the Sun is an exceedingly 
typical star in our Galaxy. Therefore, we assume that the 
frequency of stellar flares throughout the Galaxy can be 
represented by that of solar flares, i.e. the frequency of 
flares on the Sun. 

Solar flares are classified into A, B, C, M and X 
classes, by size, as identified by X-ray observations [11]. 
They are represented by a logarithmic scale such that the 
difference in intensity between classes Al and X10 is a 
factor of 100,000. The acceleration of particles over 100 
GeV is usually observed for the largest classes of solar 
flare, such as class X10 [12]. It has been reported that 
the frequency of solar flares is inversely proportional to 
the size of the flare, and the product of the frequency 
and size is constant [13]. Based on this estimation, a 
flare of size X10 occurs every 3xl0 7 seconds (once a 
year), while flares of size Ml occur every 3 days on 
average. 

(2) The energy spectrum of the accelerated protons 
for each star can be represented by an example observed 
solar flare of April 15, 2001 (the Easter event). In this 
flare the energy spectrum of the accelerated protons can 
be expressed by a power law with gamma index T = 
-3.75 (differential index) in the energy range between 
Ep = 200 MeV and 100 GeV [12]. The proton fluxes 
at the Earth for Ep >100 MeV, >1 GeV and > 10 GeV 
are 100, 1 and 10~ 3 protons/(cm 2 sec str), respectively. 
(We call these numbers f|B). 

(3) In this paper the energy spectrum of photons 
induced by stellar flares is represented by a similar 
power spectrum for the photons in the energy range 
£ 7 «17MeV to 8 GeV with spectral index T = -3.75. 
Here we have assumed the inelasticity of gamma-rays 
to be if 7 «l/12. 

A detailed Monte Carlo calculation is necessary on 
the secondary production process; however, the gross 
features of the interactions are not changed dramatically 
from this estimation. Therefore, the intensity of gamma- 
rays at Ey > 75 MeV may be described by the proton 
intensity at Ep >900 MeV — 1 GeV. It must be noted 
here that pions will be not produced by protons with 
energy less than Ep< 1 GeV. (See note in [18].) 
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(4) In the largest solar flares, the peak intensity of 
accelerated protons occurs during a period of 1200 
seconds on average. Therefore if 2.5 xlO 4 stars out of 
a total 1.6xl0 10 stars in the Galactic bulge produce 
one XlO class stellar flare once every year, high energy 
protons are constantly produced in the bulge on average 
(1200 secx2.5xl0 4 = 3xl0 7 sec = 1 year). Hence, the 
bulge is likely to be the main source for high energy 
protons of Ep > 1 GeV. 

Dividing the total number of stars in the bulge by the 
number of flare stars with intensity XlO per year, we 
obtain 1.6xl0 10 stars/2.5xl0 4 stars = 6.4 x 10 5 (we call 
this number flC), which represents the number of stars 
in the Galactic bulge that may generate a constant high 
energy proton current with energy Ep > 1 GeV. 

(5) To calculate the total high energy proton flux in- 
duced by solar flares we use the total number of protons 
accelerated from the flare of April 15, 2001 as a standard 
value. The observed flux of accelerated protons for the 
Easter flare was observed to be ~lproton/(cm 2 secstr) 
at Earth for Ep >1 GeV (from |B). The total flux at 
the Sun of the Easter flare is derived by multiplying the 
distance correction factor between the Sun and the Earth, 
i.e., 4irr 2 = 2.0xl0 23 m 2 = 2.7x 10 27 cm 2 str, where r 
corresponds to 1 AU = 1.5 x 10 11 m. In other words, if 
protons were emitted in an isotropic direction from the 
Sun, the high energy proton flux with Ep > 1 GeV (for 
which jJB~ 1) at the Sun is 2.7 xlO 27 protons/sec. This 
value corresponds to a proton current of 4.3 x 10 8 A with 
Ep >1 GeV in the flare. (We call this number jjD) 

III. Calculations and Results 

In this section, we obtain the expected flux of protons 
accelerated by stellar flares in the bulge and compare 
them with observed data. 

(1) The total high energy proton intensity produced by 
stellar flares in the bulge is determined by estimating the 
total proton current induced at the bulge using (JC. We 
assume that each star in the bulge releases a large stellar 
flare once every year, creating a strong intensity that 
continues for about 1200 seconds. We further assume 
that the flux of high energy protons generated by the 
flare can be approximated by the Easter flare. Then, 
the total flux induced by the stellar flares in the bulge 
can be represented by jjCxflD [(number of effective 
stars) x (total number of protons produced by each 
star)] or (6.4 x 10 5 ) x (2.7 x 10 27 ) = 1 .7 x 10 33 protons/sec 
for £p>l GeV. This corresponds to 2.7xl0 14 A and 
indicates that the Galactic bulge constantly generates a 
very strong DC current. 

(2) When calculating the total number of high energy 
protons in the Galaxy it is important to consider the 
Galactic magnetic field. If there was no magnetic field 
in the Galaxy, the generated protons, generated at a rate 
of 1.7 xlO 33 protons/sec, would be leaving the galaxy, 
dissipating into intergalactic space. However, the galactic 
magnetic field works to retain the generated protons in 
the Galaxy for a long time. From the measurement of 



Be 10 we estimate this time to be 20 million years [14]. 
Therefore, the total number of high energy protons in 
the Galaxy is 1.7xl0 33 x 2.0xl0 7 = 3.4xl0 40 protons. 
(We call this number jJE.) 

(3) As the produced protons travel through space, 
they undergo nuclear interactions with interstellar gas, 
as pointed out by Hayakawa, Hutchinson and Morrison. 
The probability of a nuclear interaction occurring can be 
calculated by assuming a total cross-section of 40 mb for 
protons and a gas density of one hydrogen atom/cm 3 . 
For a path of 1 kpc (= 3xl0 21 cm), the collision rate 
is 40xl0~ 27 [cm 2 ] x3xl0 21 [cm] x 1 [/cm 3 ] = 
1.2xl0~ 4 . Over 20 million years, protons with £p>l 
GeV traveling at about light speed travel a distance 
of about 6.1xl0 3 kpc, giving a product of (6.1xl0 3 ) 
x (1.2xl0~ 4 ) = 0.73. Hence, about 73 %of produced 
protons lose their energy by nuclear interactions with 
intermediate gas. (We call this number jJF). 

(4) The production rate of gamma-rays Ej>70 
MeV in the Galaxy is estimated by (JExjjF= 
3.4xl0 40 protonsx0.73 = 2.5 xlO 40 photons/sec. (We 
call this number (JG.) 

(5) For the observation of the Galactic bulge by 
gamma-ray telescopes, we estimate the geometrical ac- 
ceptance of the bulge by detectors on the Earth. The 
Galactic center is located about 8 kpc from the Earth. 
The intensity of the gamma-rays emitted at the center of 
the Bulge is reduced by the geometrical factor l/47rr 2 = 
1.3xl0- 46 crrr 2 for r = 8 kpc = 2.5 xlO 22 cm. (We call 
this geometrical reduction factor jJH.) 

(6) Since the data of the COS-B satellite is given in the 
units of (cm 2 sec str) -1 [15], we must calculate the solid 
angle of the bulge. While the solid angle of the Galactic 
center cannot be definitely determined, the bright region 
of the Galactic bulge spreads over approximately 5 
degrees in the Galactic longitude direction and 2 degrees 
in Galactic latitude. Therefore, we estimate the field 
of view of the Bulge to be 5x2 square degrees. This 
corresponds to an opening angle of 0.087 radians x 0.035 
radians or 3.0 xl0~ 3 steradians. The data of COS-B 
towards the Galactic bulge give a gamma-ray flux for the 
Galactic center of 8xl0 _4 cm" 2 sec _1 str _1 at £ 7 >70 
MeV to 5 GeV for t = 11 - 20°. Multiplying the field 
of view by this number gives an intensity of gamma- 
rays per unit area from the direction of the Bulge of 
2.4xl0~ 6 cm~ 2 sec~ 1 . (We call this number fll) 

(7) We now compare the observed data with our 
predictions. According to our prediction, the flux of 
gamma-rays from the Galactic center, obtained by 
multiplying (JG and jjH, is (2.5 xlO 40 photons/sec) 
x (1.3xlO- 46 crrT 2 ) = 3.3 x 10- 6 photons/cm 2 sec at 
£ 7 >70 MeV. (We call this number jjj.) The predicted 
number jjj is 35 % higher than the observation data 
jjl. However, given the crude estimation process, the 
coincidence of the two numbers is surprising. This shows 
that it may be possible to determine the origin of cosmic 
rays by simple calculations of this type. 

We are aware that there are numerous points of 
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possible improvement for future studies. For example, 
we need to determine if the value of 20 million years 
is reasonable for the time cosmic rays are retained in 
the Galaxy and whether the photon intensity be fairly 
approximated by the proton flux divided by a of factor 
12. 

Furthermore, we have not taken account of the con- 
tribution by primary electrons, which are expected to 
contribute to the flux of gamma-rays in the low en- 
ergy region E 1 <70MeV since electrons are simultane- 
ously accelerated with protons in solar flares. Probably 
electrons will be accelerated by stellar flares upto ~ 
100 Me V: the same Lorentz factor of protons. 

We must also take into account flares induced by 
nearby stars, which are likely to be quite strong. Addi- 
tionally, we must consider the re-acceleration process of 
seed protons produced by flares in super novae remnants. 
By these processes protons will be re-accelerated into 
higher energies ( the second step acceleration) and will 
have a harder energy spectrum. These protons should 
properly be classed as cosmic rays. 

IV. Summary and Conclusions 

(1) The origin of the diffuse gamma-rays is likely to 
be protons accelerated by stellar flares (the first step 
accleration). They are not chiefly produced by super 
novae or the interstellar medium. See Figure 2. 

(2) The energy spectrum of primary protons may be 
expressed by Ei T with the index T = -3.75 

(3) The accelerated protons by stellar flares are re- 
tained not only near the Galactic bulge but also in the 
Galactic disk and halo. They are retained in the Galaxy 
for approximately 20 million years. 

(4) When these protons collide with the interstellar 
medium, gamma-rays are produced. These gamma-rays 
are likely to be the diffuse gamma-rays observed by the 
COS-B, EGRET and Fermi satellites around £ 7 «100 
MeV [16-17]. For the production process, our model is 
based on the Hayakawa— Hutchinson— Morrison hypoth- 
esis. However contribution of pp collisions would be 
not restricted down to 70 MeV and in the lower neergy 
region, probably the contribution of electrons would be 
dominated. 

(5) The primary accelerated particles are re- 
accelerated in the interstellar shocks to the energy E^. 
The energy spectrum can be written by E\ cLE^and T=- 
2.75. 

(6) The observed flux of the diffuse gamma-rays 
with E 7 >a few GeV must be produced by these re- 
accelerated protons. However, the intensity of diffuse 
gamma-rays with around _E 7 ^70MeV may be produced 
by primary accelerated protons by stellar flares. 
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Fig. 1: A famous picture given by Lingenfelter on the 
origin of cosmic rays. In very early days, the origin 
of cosmic rays was considered as to be super nova, 
however after Fermi's idea, the origin of cosmic rays was 
considered as interstellar matter. At Kyoto conference in 
1979, Lingenfelter had pointed out that the seed particles 
of cosmic rays should be both origin named SNISM 
supernova interstellar matter 
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Fig. 2: A picture where seed particles of cosmic rays are produced. The seed particles of cosmic rays may be coming 
from stellar flares in the Galactic bulge and arm. The red cross marks represent where stellar flares occur. Those 
accelerated protons by stellar flares will be transported to the Sun through the Galactic arm. They will sometimes 
penetrate the region of super nova remnants. There must be matter turbulance and/or strong magnetic field. Then 
either by the original second Fermi acceleration mechanism or the first Fermi acceleration mechanism, these particles 
will be re-accelerated into high energies. They will have harder energy spectrum than that of original accleerated 
protons. These are called as cosmic rays. The basic picture is given by the coutesy of Cambridge University Press. 
ISBN 521 26369 7 (1985) The Cambridge Atlas of Astronmy, edited by J. Audouze and G. Israel. 



